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ABSTRACT 

We measure the zero point of the Stanek (1996) extinction map by comparing 
the observed (V — K) colors of 206 K giant stars with their intrinsic {V~K)q colors 
as derived from their H/? indices. We find that the zero point of the Stanek map 
should be changed by AAy = —0.10 ± 0.06 mag, obtaining as a bonus a three- 
fold reduction of the previous statistical error. The most direct way to test for 
systematic errors in this determination would be to conduct a parallel measurement 
based on the (V — K) colors of RR Lyraes (type ab). 

Subject Headings: Galaxy: general - Hertzsprung-Russell diagram 
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1. Introduction 



Baade's Window, {£,b) ~ (1°, —4°), has been an important laboratory for the 
study of bulge populations. The key features of Baade's Window that made it 
such a focus of early work were its relatively low extinction (Ay ~ 1.5) and the 
presence of NGC 6522 which provided an opportunity to measure that extinction. 
The value of this window increased significantly when Stanek (1996) constructed 
a detailed extinction map with 30" resolution by applying the method of Wozniak 
& Stanek (1996) to observations by OGLE (Szymahski & Udalski 1993; Udalski 
et al. 1993) over a 40' square field. It is now possible to dercdden stellar samples 
on almost a star-by-star basis in this field. Stanek (1996) estimates the error in 
differential extinction to be ~ 0.1 mag in Ay, but notes that the errors in the 
absolute extinction are dominated by the zero-point error, 0.20 mag. For many 
applications, such as the interpretation of color-magnitude diagrams of bulge field 
stars and of the cluster NGC 6522 or the measurement of distances using RR 
Lyraes or other tracers, the determination of the zero point is crucial. 

Here we argue that the best way to estimate the zero point for the Ay map is 
to measure AE{V — K)i for an ensemble of stars i = l...n, defined by 

/^E{V - K)i = {V- K)i -{V- K)o,i - (1 - a)4;r^ (1.1) 

where {V — K)i is the observed color of the star, [V — -ft')o,i is its predicted unrcd- 
dened color, ^Stanek jg ^j^g visual extinction at the position of the star in the Stanek 
(1996) map, and a is the ratio of the extinction in the K and V bands, assumed 
to be (Rieke & Lebofsky 1985), 

Ak , . 

a =-^ = 0.11. (1.2) 

The correction to the zero point of the Stanek (1996) map is then given by 

= <^y-^)) , (1.3) 

where (A£'(y — K)) is a suitably weighted average of equation (1.1) over the 
sample. We then apply this method to the data of Terndrup, Sadler, & Rich 
(1995, TSR) and find 

AAy = -0.10 ±0.06. (1.4) 

That is, the Ay of stars in Baade's Window are on average 0.10 mag Zow;er than 
previously believed. 
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2. The Stanek Extinction Map 

Stanck (1996) measured the mean position of the red giant clump on the color- 
magnitude diagram {{V), {V — I)) as a function of position on the sky {x,y). He 
then estimated the differential total and differential selective extinctions, 

Av{x, y) = {V){x, y) + Cy, E{V - I){x, y) ^ {V - I){x, y) + Cy-i, (2.1) 

with the constants, Cy and C^/-/) being undetermined. He found empirically a 
very strong correlation between Ay and E{y — I): 

Av{x, y) ~ 2.49 E{V - I){x, y) + C. (2.2) 

On physical grounds, C = 0. This left one constant, either Cy or Cy_/, to be 
determined. 

Stanek (1996) then determined the overall zero point using the measurement 
E{y — I) = 0.59 ± 0.08 made by TSR in a sub-region of Baade's Window, Blanco 
region A (Blanco, McCarthy & Blanco 1984). That is, he set the zero point so that 
the average of E(y — I) over this region of his map reproduced the TSR value. He 
then used equation (2.2) (with C = 0) to establish the zero point of the Ay map; 
i.e., he set the mean extinction of region A to be {Ay) = 0.59 x 2.49 = 1.47. The 
formal error in this determination is therefore 0.08 x 2.49 0.20. 

3. Previous Approaches 

With the exception of TSR, all previous determinations of the extinction 
toward Baade's Window have been made by measuring the selective extinction 
E[B — V), and then multiplying by an assumed ratio of total to selective extinc- 
tion Ry = Ay/E{B - V) (Arp 1965; van den Bergh 1971; Walker & Mack 1986; 
Terndrup & Walker 1994). There are several major disadvantages to this approach. 
First, Ry ~ 3 is rather large, and the statistical error in E{B — V) (usually es- 
timated to be > 0.03) is multiplied by this factor when estimating the error in 
Ay. Second, Ry varies along different lines of sight, so for any particular hne of 
sight for which it is not actually measured (e.g., Baade's Window) the precision 
of the estimate is no better than 7%. Hence, the statistical error alone for Ay is 
more than 0.12 mag. Finally, there are systematic errors arising from uncertainties 
in the intrinsic {B — V) colors of stars used to estimate E{B — V). While the 
intrinsic color of extremely hot stars (in the Raleigh- Jeans limit) is known from 
fundamental physics, there are no such stars lying beyond the dust column in this 
direction. Hence, one must use cooler stars whose {B — V) colors are sensitive 
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functions of temperature, metallicity, and perhaps other factors. The standard 
approach is to find local analogs of the program stars and directly measure their 
colors, but systematic errors may arise from any unrecognized differences between 
these two groups of stars. As always, it is difficult to determine the size of the 
systematic errors, but one can gain a sense of their magnitude by comparing the 
E{B — V)q = 0.45 ± 0.04 derived by van den Bcrgh (1971) from three different 
methods based on cool stars (K and M giants) with the E{B — V)o = 0.60 ± 0.03 
derived by Walker & Mack (1986) using relatively hot stars (RRab Lyraes). Here 
the subscript "0" means "reduced to zero color" using the prescription adopted by 
TSR from Dean, Warren, & Cousins (1978). 

TSR pioneered a radically different approach, although they did not call explicit 
attention to this fact. They measured 

£;(!/- X) = 1.23 ± 0.08, (3.1) 

for the Blanco A region by comparing the if/3A4861 (Faber et al. 1985) index as 
a function of calculated {V — K) color to that observed for bright K giants in the 
solar neighborhood. They then inferred (but did not explicitly write down), 

Av ^ ^ 1.38 ± 0.09. (3.2) 

1 — a 

While this approach is formally identical to the previous one (measurement of 
a selective extinction and conversion to a total extinction), it is potentially more 
accurate than using E{B—V) because the extrapolation to total visual extinction is 
small (a factor 1.12 vs. 3), and therefore the error in Ay is only slightly bigger than 
the error in E(y — K). TSR then used this measurement to infer E{B — V)q = 
0.47 ± 0.04 and E{V - I) = 0.59 ± 0.08. To obtain these quantities, they had 
to employ estimates of the ratios of total to selective extinctions which accounts 
for the proportionately larger error bars compared to those in equations (3.1) and 
(3.2). In particular, TSR assumed 

Ay = 2.33E{V - I), (3.3) 

where the coefficient is considerably lower than Stanek's (1996) empirical value for 
Baade's Window [eq. (2.2)]. By fixing the zero point according to TSR's E{V — I) 
rather than TSR's more directly determined Ay, Stanek (1996) therefore overes- 
timated the visual extinction by (2.49 — 2.33) x 0.59 = 0.09 mag, and also overes- 
timated the uncertainty. Prom this simple argument, we therefore derive a naive 
correction to the Stanek (1996) extinction map, 

AAy = -0.09 ± 0.09 (naive). (3.4) 
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4. New Determination 



However, rather than simply adopting the correction given by equation (3.4), 
we prefer for several reasons to make a complete redetermination of the zero point 
of the Stanek (1996) extinction map. First, Stanek (1996) compared the TSR 
extinction for the whole of region A with the mean extinction of the subregion of 
region A that is covered by his map. About 35% of the TSR stars he within 2' of 
NGC 6522 and so are excluded from the map. One should therefore compare the 
mean extinction of the TSR stars lying within the map with the mean extinction 
predicted by the map for their positions, otherwise systematic trends of extinction 
with position could affect the result. Second, TSR measured and reported the mean 
extinctions in region A and in regions B and C separately. They did so because 
these areas have different average extinctions (Blanco et al. 1984), so it would 
degrade the information content of their results to combine the two. However, for 
purposes of measuring the offset to the Stanek (1996) map, the fact that different 
regions have different extinctions is irrelevant. The only concern is to measure the 
difference in the observed extinction at each point from that predicted by the map. 
Including all three regions approximately triples the sample and correspondingly 
reduces the statistical errors. Third, TSR used the relatively crudely determined 
}if3/(y—K)Q relation of Faber et al. (1985). A much more sophisticated polynomial 
relation is available from Gorgas et al. (1993, Table 6, eq. 5): 

H/3 =8.2261 - 5.9295(y - K)o + 0.52968(logc/) - 0.048352(log5)2 

- 0.23695(log g) {V - K)q + 1.8169(y - K)l - 0.19721(y - K)l, ^^'"^^ 

where we have ignored the terms in [Fe/H] since they are small and since ([Fe/H]) ~ 
for the sample in any case. 

To carry out our analysis, we first restrict the sample to the 209 stars with 

2.1 < (y - X), < 3.0, (4.2) 

where 

{V - K), = {y-K)-{l- a)(^5^^^«k ^ ^^^^^^^ ^4 3^ 

is our final best-fit estimate of {V — K)q as determined from our final best-fit offset 
to the Stanek (1996) map [eq. (1.4)]. The upper limit is chosen to exclude M 
giants which have TiO bands that influence the Hj3 index. Gorgas et al. (1993) 
also excluded these stars when they derived equation (4.1). We note, however, 
that if we extend the sample to the 239 with (y — X)* < 3.8 (the range of the 
Gorgas et al. 1993 data), the results change by less than 0.01 mag. The lower limit 
in equation (4.2) is justified below. 
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To obtain (V — K) from the V and / measurements, we slightly modify the 
procedure of TSR. TSR converted from (V — I) to (V — K) colors based on the 
extremely tight quadratic color-color relation obtained for a patch of Baade's Win- 
dow by Tiede, Frogel, & Terndrup (1995). We modify this procedure by using the 
Stanek (1996) differential map to transform each {V — I) from its observed value 
to the value it would have if it lay in the Tiede et al. (1995) region and then use 
the same relative extinction to transform the inferred {V — K) back to the star's 
actual position. In practice, the individual corrections generated by this procedure 
are small, typically < 0.01 mag, because the color-color track is almost parallel to 
the reddening vector. The net effect on the final result is <^ 0.01 mag. 

We estimate the surface gravity g from the inverse square of the star's effective 
radius in K band, 

s 

log A = log + o.4(X - aAv - Mk,q - //) , (4.4) 

where S oc [exp{hc/ XxkT) — is the black body surface brightness at K band 
{Xk = 2.2^m) for an assumed temperature T = 8520-2230(l/-ir)o+267(\/-ir)^ 
obtained by fitting the values given in Table 4 of Ridgway et al. (1980). Here 
jj, = 14.5 is the adopted distance modulus to the Galactic center and Mx,q = 3.3. 

For each trial value of AA\/, we deredden the "observed" {V —K) color (inferred 
from V — I) and K magnitude of each star, use these to estimate its temperature 
and surface gravity, and finally predict H/? from equation (4.1). We form x^(AAy) 
from the difference of the observed H/3 and these predictions divided by the errors 
as reported by TSR. The best-fit AAy and its errors are determined from this 
function. 

We find that in the adopted interval, 2.1 < {V — K)^ < 3.0, the predicted 
and observed H/9 indices (weighted by the observational errors) are in good overall 
agreement with = 179 for 208 degrees of freedom. However, the observations 
deviate markedly from the predictions for {V — K)^ < 2.1 (see Fig. 1, below). 
These bluer stars are mostly G giants and subgiants. The original sample from 
which TSR drew their stars was selected primarily for proper motion studies and 
hence was composed of preferentially brighter stars. The intrinsically fainter G 
stars are therefore likely to be foreground disk stars. From equations (4.1) and 
(4.4), we find that if one of these star lies A/i in the foreground, H/? will be 
overestimated by ~ 0.1 A/x. Thus, this selection effect can explain at least some of 
the observed deviation. In any event, these stars clearly differ from the bulk of the 
K giants (which dominated the fit by Gorgas et al. 1993) and we therefore restrict 
the sample to (V — X)* > 2.1. 
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We now consider various systematic effects. First, there arc several assumptions 
that affect the results through the estimate of the surface gravity. For example, we 
find that if the adopted distance to the Galactic center is increased by A//, then 
the surface gravities are decreased, implying that AAy is also decreased (becomes 
more negative) by — O.llAyU. Similarly, if the mean mass of the stars is larger than 
one solar mass (the value that we implicitly assumed in writing cq. (4.4)) then 
AAy is increased by O.lOAlnM. If the effective radius at which surface gravities 
are measured differs from that of the K band photosphere, AAy is changed by 
— 0.20Alni?. If the adopted temperatures are on average different from the true 
ones, A Ay is changed by 0.19 A InT. Plausible adjustments for each of these factors 
are therefore likely to affect AAy by ^ 0.01 mag. 

Gorgas et al. (1993) report a scatter in their fit to the H/3 index of 0.28. We 
estimate that this scatter causes an uncertainty in the zero point of equation (4.1) 
of ~ 0.025, which translates directly into an uncertainty in AAy of 0.03 mag. 

Next, we break the sample into two subsamples, one of stars in region A (where 
TSR found Ay = 1.38) and one of stars in regions B and C (where they found Ay — 
1.55). We obtain separate fits of AAy = -0.07 ± 0.09 and AAy = -0.13 ± 0.06, 
which are consistent at the 1 a level. 

Finally, we divide the sample into two subsamples, according to their angular 
distance 9 from the center of NGG 6522. The inner group (2.'0 < 9 < 2.'5) has 

49 stars and the outer group {6 > 2.'5) has 160 stars. We find values of AAy = 
—0.28 ±0.12 and AAy = —0.06 ±0.06, respectively. That is, they are inconsistent 
at the 1.6(7 level. This is somewhat worrisome because it may indicate that the 
inner group is affected in some way by proximity to the cluster. We investigate 
the following possible effects. First, Stanek (1996) reports that the extinction 
values close to NGG 6522 and NGG 6528 are systematically lower than surrounding 
regions, leading him to believe that he may have underestimated it (due to some 
unspecified form of contamination). This was the primary reason for excluding the 
region 6 < 2' from the map. It is possible that this effect, if real, extends beyond 
2'. However, the sign of the effect is wrong to explain the difference between the 
two subsamples. It is possible that contamination by cluster stars generates some 
other effect that has the correct sign. We perform the following tests to search for 
cluster contaminants. First we search for an excess of stars in the underlying sample 
(which reaches as close as ~ 1' from the cluster center) with radial velocities that 
are consistent with the cluster velocity, ~ — 25kms^^ (Rich 1990; Smith, Hesser & 
Shawl 1976), both in the sample as a whole and as a function of 9. If contamination 
by the cluster were significant, one would expect an excess, especially at small radii. 
None is detected. Next we conduct a test that is sensitive to even lower levels of 
contamination: we plot the proper motions of all stars having radial velocities 
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consistent with cluster membership. If even a small subset of these are in the 
cluster, the proper-motion diagram should show a clump. We detect a common 
proper motion clump of eight stars. We will report elsewhere on this measurement 
of the proper motion of NGC 6522. For present purposes, we note that three of 
the eight stars are part of our sample, all three being in the inner annulus. We 
exclude these, leaving a sample of 206 stars. The estimates for AAy in the inner 
and outer annuli are then —0.25 ± 0.12 and AAy = —0.06 ± 0.06, respectively, 
a 1.4(7 difference. Since contamination by the cluster has been eliminated, this 
difference should be regarded as a normal statistical fluctuation, and we therefore 
include the entire remaining sample of 206 stars and flnd 

AAv = -0.10 ± 0.05 (internal error). (4.5) 

We then add in quadrature the external calibration error of 0.03 mag intrinsic to 
the Gorgas et al. relation (4.1) to obtain equation (1.4). Figure 1, shows the mean 
predicted and observed H/3 indices (weighted by the observational errors) for the 
overall best fit, binned by {V — K)^. We find — 176 for 205 degrees of freedom. 

For completeness we note that had we adopted the relation of Faber et al. 
(1985) in place of equation (4.1) from Gorgas et al. (1993), AAy would decrease 
from -0.10 to -0.11. 

5. Future Tests 

Of course, the most important potential source of systematic errors is not 
probed by the tests of the previous section: the possibility that the stars in Baade's 
Window differ systematically in some unknown way from the local stars upon 
which equation (4.1) is based. Ultimately, the only way to test for this effect is 
to make an independent determination of E{V — K) on a. substantially different 
set of stars. The obvious choice for this test is RR Lyraes (type ab). First, the 
comparison between RRab's and cool giants showed the largest discrepancy of all 
determinations based on E{B ~ V), so it is important to see if this discrepancy 
persists for E{V — K). Second, because of the relatively narrow range of RRab 
colors and the accuracy of the Stanek (1996) map, we estimate that each star 
should provide a statistically independent estimate of the zero point accurate to 
0.15 magnitudes. Since there are more than 50 such stars in the region, the limit 
for this method is set by the size of the calibrating sample (17) observed by Jones 
et al. (1992). We estimate this limiting uncertainty to be only ±0.03 mag. Work 
is in progress to apply this method. 
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95-28227 from the NSF. 
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FIGURE CAPTIONS 



1) Observed H/5 index {solid squares) compared to the values predicted {open 
circles) from the relation of Gorgas et al. (1993). The prediction is mainly a 
function of {V — K)^, the best-fit estimate of (V" — K)q, using the observed 
{V — K) color, the Stanek (1996) extinction for the position of each star, 
and the overall offset given by eq. (1.4). The points are binned by ~ 23 
stars and the error bars indicate the standard error of the mean based on the 
individual errors reported by TSR. The stars {V — K)^ > 3.0 (right dashed 
line) were excluded from the fit because they are M giants which are affected 
by TiO bands and because Gorgas et al. (1993) excluded such stars when 
they determined eq. (4.1). The stars {V — K)^ < 2.1 (left dashed line) were 
excluded from the fit because they are contaminated with foreground stars 
and hence fall systematically below the predicted values (see text). 
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